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SUMMARY 

INTRODUCTION 

Computer  programs  often  require  a means  of  expressing  functional 
relationships  in  which  no  exact  mathematical  expressions  are  known.  In 
this  situation  tabular  data  inputs  can  be  used  for  these  functional 
relationships.  Numerous  techniques  are  available  for  interrogating 
table  data  arrays,  from  the  simplicity  of  linear  technique  to  the  complexity 
of  higher  order  polynominal  techniques.  The  purpose  of  this  report  is 
to  describe  a complete  computer  routine  package,  TLOOK/ TREAD,  for  input 
and  interpolation  of  tabular  data  which  is  accurate  and  has  broad  application 
to  both  the  vehicle  design  and  technology  programs  at  the  Naval  Air 
Development  Center  as  well  as  other  programs  which  require  table  look-up 
information.  The  computer  code  can  determine  the  value  of  a dependent 
variable  in  terms  of  up  to  three  independent  variables  using  a cubic 
spline  interpolation.  The  inputs  to  this  code  require  an  effective 
means  of  locating  and  correcting  input  errors.  Another  report  describing 
a companion  computer  code,  GPPR  (General  Purpose  Plotting  Routine), 
which  produces  CALCOMP  plots  of  tabular  data  processed  by  TLOOK/TREAD 
will  shortly  be  published  as  reference  (a). 

SUMMARY  OF  RESULTS 

A general  purpose  table  look-up  cubic  spline  interpolation  technique 
has  been  programmed  for  the  CDC  6600  computer.  The  utilization  of  the 
resulting  package  of  routines  (called  TREAD/TLOOK  and  SPLNQ1)  is  illustrated 
with  a sample  problem.  A comparison  between  actual  and  interpolated 
values  of  dependent  variables  in  the  sample  problem  has  shown  that  this 
interpolation  technique  is  a fast,  accurate  and  a valuable  computer 
tool.  The  programmed  cubic  spline  interpolation  technique  used  by 
TREAD/TLOOK  is  greatly  improved  over  that  shown  in  reference  (b).  The 
major  improvements  concern  computational  speed,  computer  storage  size 
and  flexibility.  The  input  format  to  TREAD/TLOOK  is  arranged  to  accept 
both  symmetrical  and  nonsymmetrical  tabular  data  arrays. 
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DISCUSSION 


BACKGROUND 

Numerous  codes  in  current  use  at  NAVAIRDEVCEN  require  interpolation 
schemes  as  a means  of  expressing  functional  relationships  between  dependent 
and  independent  variables.  Linear  and  higher  order  polynominal  curve 
fits  represent  two  approaches  for  this  purpose.  Linear  interpolation 
schemes  have  the  advantage  of  simplicity  and  fast  computational  speed. 

A search  procedure  is  used  to  determine  the  correct  data  point  interval 
where  the  interpolation  is  performed.  However,  unless  a large  number  of 
data  points  is  input,  interpolation  accuracy  may  be  lost.  In  addition, 
the  first  derivative  of  the  linear  interpolated  function  is  discontinuous 
at  each  input  data  point.  Higher  order  polynominal  interpolation  schemes 
exhibit  some  of  the  same  inaccuracies  as  linear  interpolation  schemes. 

The  function  obtained  from  the  interpolation  may  or  may  not  coincide 
with  each  input  data  point.  In  addition,  curve  fits  tend  to  oscillate 
about  the  data  points  resulting  in  numerous  inflection  points.  However, 
this  approach  has  the  advantage  of  eliminating  the  search  procedure 
necessary  for  linear  interpolation  schemes.  The  approach  suggested  in 
this  report  circumvents  some  of  these  disadvantages. 

PROGRAM  PACKAGE  DESCRIPTION 

The  table  look-up  program  package  in  this  report  consists  of  the 
TREAD/TLOOK  routine  and  function  SPLNQ1.  TREAD/TLOOK  is  the  routine  by 
which  tabular  data  is  read  as  input  and  stored  for  the  spline  fit  interpolation 
by  the  SPLNQ1  routine.  A user's  guide  in  Appendix  A illustrates  the  use 
of  the  table  look-up  package.  A general  description  and  operation  of 
the  TREAD/TLOOK  and  SPLNQ1  routines  follows. 

Description  of  Subroutine  TREAD/TLOOK 

Subroutine  TREAD  (with  entry  TLOOK)  is  programmed  to  interpolate  or 
extrapolate  a given  set  of  data  with  one,  two,  or  three  independent 
variables.  The  actual  interpolation  or  extrapolation  is  performed  by 
the  function  SPLNQ1  and  will  be  discussed  later.  The  TREAD  portion  of 
the  subroutine  reads  input  data  and  assembles  the  data  in  the  format 
required  for  function  SPLNQ1.  The  TREAD/TLOOK  subroutine  has  the  capability 
to  accept  up  to  30  sets  of  data,  each  set  representing  a table  look-up 
function  of  up  to  three  independent  variables.  The  user  initiates  a 
TREAD  operation  to  input  table  data.  As  the  input  data  is  processed  and 
stored  in  a single  dimension  array  (called  TDATA)  the  starting  array 
location  and  table  reference  number  of  each  table  is  also  stored.  The 
table  reference  number  is  a 1 to  4 digit  number  used  to  identify  each 
table.  To  use  a particular  table  look-up  function  the  user  merely 
specifies  the  appropriate  table  reference  number  and  the  values  of  the 
independent  variables  in  a call  to  the  TLOOK  routine.  When  a TLOOK 
operation  is  initiated  a search  of  the  table  reference  numbers  is  performed 
to  determine  the  TDATA  array  location  of  the  table  data.  Once  the 
initial  TDATA  array  location  is  determined,  the  function  routine  SPLNQ1 
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Is  used  repeatedly  to  interpolate  and/or  extrapolate  the  table  data. 

The  first  30  array  locations  in  the  TDATA  array  are  reserved  for  values 
of  table  reference  numbers  and  the  next  30  array  locations  are  reserved 
for  values  indicating  the  starting  location  of  each  table  data  set.  The 
TDATA  storage  locations  starting  at  61  are  reserved  for  the  actual  table 
data  in  groups  of  3n+3  storage  locations  as  illustrated  in  the  example 
which  follows. 

Example 

In  Figure  1 a three  independent  variables  table  look-up  function  is 
shown.  The  first,  second  and  third  independent  variables  are  respectively 
X,  Y,  and  Z.  In  this  case  Z represents  a plane  and  X and  Y represent 
independent  variables  lying  on  any  plane,  Z = constant.  The  dependent 
variable  is  FXYZ.  The  object  of  the  interpolation  and/or  extrapolation 
is  to  determine  FXYZ  from  any  set  of  input  values  Xi,  Yi,  and  Zi.  In 
this  example  the  TDATA  array  would  contain  the  following  3n+3  groups: 

Z data 
Y data 
X data 
FXYZ  data 
X data 
FXYZ  data 


Y data 
X data 
FXYZ  data 
X data 
FXYZ  data 


Y data 
X data 
FXYZ  data 
X data 

The  final  value  of  a three  independent  variable  interpolation  is 
obtained  after  three  steps.  The  first  step  is  to  interpolate  or  extrapolate 
(in  each  constant  Z plane  and  along  each  constant  Y line)  in  terms  of 
the  first  independent  variable  X.  This  step  yields  FXYZ  for  the  given 
value  of  Xj.  The  second  step  is  to  interpolate  or  extrapolate  (in  each 
constant  Z plane)  the  values  of  FXYZ  (for  the  given  values  of  X ^ ) in 
terms  of  the  second  independent  variable  Yi.  This  step  yields  FXYZ  for 
the  given  value  of  Xi  and  for  the  given  value  of  Yi.  The  final  step  is 
merely  a single  interpolation  or  extrapolation  of  3n+3  sets  of  data 
stored  in  a single  dimension  array  in  terms  of  the  third  independent 
variable  Zi. 
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In  each  set  of  3n+3  storage  locations  the  following  information  is 
stored: 

1.  the  number  of  data  points,  n 

2.  the  ascending  order  values  of  the  independent  variable  at  each 
data  point 

3.  the  corresponding  values  of  the  dependent  variable  at  each  data 
point 

4.  a value  indicating  which  two  adjacent  data  points  were  used  in 
the  last  interpolation 

5.  a value  indicating  either  the  last  calculated  second  derivatives 
can  be  used  again  or  new  second  derivatives  must  be  calculated 

6.  n storage  locations  reserved  for  the  computed  second  derivatives 
at  each  data  point 

A FORTRAN  listing  and  brief  description  of  the  calculation  sequence 
in  TREAD/TLOOK  is  found  in  Appendix  D. 

SPLINE  Interpolation  Technique 

An  efficient  method  of  interpolating  is  to  plot  data  points  and  draw 
a smooth  curve  through  these  using  a draftsman's  spline.  The  interpolation 
is  then  performed  by  reading  values  of  the  dependent  variable  at  the 
appropriate  point  along  the  smooth  drawn  curve.  This  same  spline  inter- 
polation technique  is  readily  programmed  on  a computer.  Basically, 
sections  of  a third  degree  polynominal  are  used  to  connect  each  pair  of 
adjacent  data  points  in  a manner  such  that  the  first  and  second  derivatives 
are  continuous  at  each  input  data  point.  The  resulting  spline  curve  is 
not  only  smooth  and  continuous,  but  also  passes  through  each  input  data 
point.  An  energy  analogy  of  the  spline  curve  is  as  follows:  An  infinitely 

thin  draftsman's  spline  is  contrained  to  pass  through  each  data  point  in 
a manner  such  that  the  strain  energy  of  the  spline  is  a minimum.  The 
computer  solution  of  the  spline  fit  equation  involves  solving  a set  of  n 
equations  with  n unknowns  using  a matrix  elimination  method,  where  n is 
the  number  of  data  points.  The  n unknowns  determined  in  this  method  are 
the  second  derivatives  of  the  spline  curve  at  each  respective  data 
point.  After  the  second  derivatives  are  found  a search  is  performed  to 
determine  between  which  two  adjacent  data  points  the  given  independent 
variable  lies.  The  interpolated  value  of  the  dependent  variable  for  the 
given  independent  variable  value  is  then  computed  using  the  following 
equations: 

y = Cj  * bk3  + C2  * (X  - Xk)3  + C3  * bk  + C4  * (X  - Xk) 
where:  y = interpolated  value  or  dependent  variable 

X = given  value  of  independent  variable 

k = subscript  of  data  point  with  value  less  than  X (determined 
from  search) 
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bk  - *k+l  - X 

dk  = Xk+i  - Xk 

Ci  = yk  / 6 dk 

C2  * y'k+i/6dk 

c3  = yk/dk  ~ ykdk/6 

C4  = yk+i/dk  - yk+i  dk/6 

A more  complete  derivation  of  the  cubic  spline  fit  equations  is 
found  in  reference  (b).  The  calculations  for  the  interpolated  value  are 
lengthy  and  complex  and  require  computations  of  second  derivatives.  For 
this  reason  considerable  computation  time  is  saved  if  the  second  derivatives 
are  saved  for  the  next  interpolation.  These  second  derivatives  may  be 
used  again  only  if  the  values  of  the  input  data  points  do  not  change. 
Provisions  for  saving  these  second  derivatives  are  included  in  SPLNQ1 
and  represents  an  improvement  over  the  programming  technique  of  reference  (b). 

The  spline  fit  curve  can  give  very  accurate  values  for  first  and 
second  derivatives,  providing  the  spline  points  are  accurately  given. 

The  derivatives  are  very  sensitive  to  small  errors  in  coordinates  when 
points  are  close  together.  The  best  procedure  is  to  use  as  few  points 
as  is  reasonable.  The  point  spacing  should  be  proportional  to  the 
curvature  and  to  the  slope  of  the  fitted  function.  With  a very  sharp 
bend  or  with  a steep  slope  near  vertical  within  a few  degrees,  spline 
points  should  be  quite  close  and  very  accurate. 

A FORTRAN  listing  and  brief  description  of  the  calculation  sequence 
in  SPLNQ1  is  found  in  Appendix  E. 
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recommendations 

It  is  recommended  that  TREAD/TLOOK  be  used  where  applicable  in  all 
future  computer  code  developments.  It  Is  further  recommended  that  the 
TREAD/TLOOK  routine  package  be  gradually  introduced  into  existing  computer 
codes.  This  will  permit  compatible  table  data  inputs  between  the  various 
operating  programs  as  required. 


REFERENCES 

(a)  Caddy,  M.  J.,  NADC  Report  No.  NADC- 76367- 30,  GPPR  - A Multipurpose 
Computer  Code  for  Data  Plotting,  unpublished 

(b)  Pennington,  R.  H. , 1970,  "Introductory  Computer  Methods  and 
Numerical  Analysis",  the  MacMillan  Company,  London 


NADC-76366-30 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


10 


NA DC-76366-30 


APPENDIX  A 


USER' S GUIDE 


NA DC-76366-30 


it 

E. 


APPENDIX  A 
USER*  S GUIDE 

Subroutine  TREAD  (with  entry  TLOOK)  is  capable  of  accepting  up  to  30 
table  look  up  functions.  The  data  for  these  table  look  up  functions  is  read 
from  input  data  cards  into  the  subroutine.  The  input  data  read-in  mode  is 
initiated  by  the  user's  program  with  a FORTRAN  statement  as  follows: 

CALL  TREAD  (II,  X,  Y,  Z,  FXYZ) 

Only  the  parameter,  II,  need  be  defined  in  the  input  data  read-in  mode.  If 
II  ■ 1 table  input  data  will  be  output.  If  II  - 0 table  input  data  will  not 
be  output.  The  table  input  data  card  deck  consist  of  groups  of  cards,  each 
group  representing  the  inputs  to  a single  table  look-up  function.  Each  group 
of  cards  representing  a single  table  look-up  function  is  composed  of  the 
following  cards: 


Card 

Item 

Variable 

Name 

Format 

Definition 

1 

1 

ITABNO 

I. 

(columns  2-5) 

table  reference  number  of 
table  look-up  function.  If 
ITABNO  ■ 0 the  table  data 
read-in  mode  is  terminated. 

Card 

Item 

Variable 

Format 

Definition 

Name 

1 

2 

ALPHA/ 

NUMERIC 

(columns  6 to  75) 

table  title  identification 

Card 

Item 

Variable 

Name 

Format 

Definition 

1 

3 

NPC 

15 

(columns  76-80) 

table  input  card  format  switch 

If  NPC  • 0 the  card  input 
format  is  (A4,  13,  7X,  8F7.0) 
If  NPC  ■ 1 the  card  input 
format  is  (A4,  13,  3X,  7F10.0) 
(see  card  2 item  3 and  on) 


Card 

Item 

Variable 

Format 

Definition 

2 

1 

ID 

A4  (columns  1-4) 

ID  is  a 4 character  identifier 
used  to  identify  the  third 
independent  variable.  If 
table  look-up  is  2 degree  or 
less  use  a dummy  identifier. 

Card 

Item 

Variable 

Format 

Definition 

2 

2 

N 

13  (columns  5-7) 

N is  the  number  of  values  to 
be  read  in  on  the  remainder  of 
this  card  and  other  cards  if 
necessary.  (N  must  be  less 
than  100) 

A -2 
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Card  Itema  Variable  Format  Definition 

2 3 and  on  A 8F7.0 

(columns  15-21, 

22-28,  etc.)  or 
7F10.0  (Columns 
11-20,  21-30,  etc.) 
see  item  3 on  card  1 

Note:  If  the  value  of  N is  greater  than  7 (when  NPC  » 0)  or  greater  than  8 
(when  NPC  ■ 1)  the  remaining  values  of  the  third  Independent  variable  (not 
Included  on  the  previous  card)  are  input  using  either  the  format  (14X,  8F7.0) 
or  10X,  7F10.0)  according  to  the  value  of  NPC,  item  3 card  1. 


Item  Variable  Format  Definition 

(All  remaining  cards  have  the  same  format  as  card  2.  The  item  which 
distinguishes  the  card  types  is  the  value  of  the  ID  variable.  The 
4 character  values  of  each  ID  variable  on  the  first  four  cards  must 
not  be  identical.  The  4 characters  of  each  ID  variable  on  the  four 
cards  (after  the  title  card)  are  user  selected.  The  card  order  of 
each  ID  variable  is  significant.  The  first  four  cards  with  respective 
ID  variable  are  as  follows: 


Card*  Definition 

2 third  independent  variable, 

ID  and  values 

3 second  independent  variable, 

ID  and  values 

4 first  independent  variable, 

ID  and  values 

5 dependent  variable,  ID  and 

values 

The  rasMlnlng  input  cards  use  these  same  ID  values  as  input  above.  On 
cards  4 and  5 are  the  dependent  and  first  Independent  variable  values  along 
the  line  given  by  the  first  value  of  the  second  independent  variable  and  in 
the  plane  of  the  first  value  of  the  third  Independent  variable. 

Cards  with  the  ssm  respective  ID  value  as  cards  4 and  5 are  repeated  for 
different  valuee  of  second  Independent  variable  until  all  second  Independent 
variables  have  been  exhausted.  The  next  card  has  an  ID  value  corresponding 
to  the  second  Independent  variable  and  new  values  of  that  variable  for  the 
plane  of  the  second  value  of  third  Independent  variable.  The  values  of  the  first 
Independent  variable  need  not  be  repeated  if  they  are  the  same  along  each 
second  Independent  equals  constant  line.  In  each  Instance  where  the  values 
are  changed  a new  card  is  required.  Each  table  is  ended  by  the  ID  variable 


♦(after  title  card) 


Card 

3 and 
following 


These  are  the  values  of 
third  independent  variable 
in  ascending  order. 
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with  characters  set  to  EOT  (end  of  table).  An  example  of  a synnetrical  table 
and  asymmetrical  table  Is  illustrated  in  the  sample  problem,  Appendix  B. 
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APPENDIX  B 
SAMPLE  PROBLEMS 


In  Appendix  B the  use  of  subroutine  TREAD/TLOOK  will  be  illustrated 
through  a sample  problem  with  two  table  data  sets. 

Table  Input  Data  Preparation 


Part  (a)  Unsymmetrical  Example 

A french  curve  was  used  to  arbitrarily  draw  the  three  curves  shown 
in  Figure  B-l.  Data  points  were  obtained  from  these  curves  at  the 
points  indicated  by  the  symbols  on  Figure  B-l.  These  data  were  then 
input  into  the  TREAD/TLOOK  subroutine.  The  input  data  cards  for  this 
table  are  shown  in  Table  B-I,  cards  1 to  10.  On  card  1,  the  table 
reference  number  and  table  identification  are  shown.  On  card  2,  the 
third  independent  variable  identifier,  number  of  values  and  values  of 
the  variable  are  shown.  In  this  table,  the  third  Independent  variable  is 
not  used  and  therefore  Z is  a dummy  identifier.  On  card  3 the  information 
for  the  second  independent  variable  is  shown.  As  indicated  on  card  3, 

ARG2  (the  second  independent  variable)  has  3 values:  1.,  2.,  and  3., 
respectively.  On  cards  4 to  9,  the  values  of  the  first  independent 
variable  (ARG1),  and  respective  values  of  the  dependent  variable  (FUN) 
are  shown.  For  example,  for  the  curve  ARG2  = 1.  shown  in  Figure  B-l, 
five  values  of  ARG1  and  FUN  are  input  on  cards  4 and  5,  respectively. 

Other  values  for  ARG2  = 2.  and  ARG3  are  input  on  lines  6 through  8, 
respectively. 

Part  (b)  Symmetrical  Example 

On  cards  11  to  32  input  cards  to  a second  table  are  shown.  This 
example  could  represent  an  aircraft  drag  polar  which  is  a function  of 
Mach  number  and  sweep  angle.  On  card  11  the  number  1 is  shown  in  column 
80  of  the  input  card  indicating  the  optional  input  data  format  is  used. 

For  this  table  SWEP,  MACH,  and  CL  are  the  third,  second,  and  first 
independent  variables,  respectively,  and  CD  is  the  dependent  variable. 

On  card  21,  CL  values  are  input  and  retained  as  the  first  independent 
variable  values  for  all  of  the  remaining  CD  values.  The  CL  values  are 
retained  because  the  CL  values  are  not  redefined  in  the  remaining  data 
cards.  This  illustrated  the  symmetrical  feature  of  the  input  scheme. 

The  MACH  values  for  the  second  value  of  the  third  independent  variable 
are  shown  on  card  28.  As  previously  described  in  the  user's  guide, 
Appendix  A,  the  second  independent  variable  values  (MACH  in  this  example) 
must  be  repeated  in  both  symmetrical  and  unsymmetrical  table  inputs. 

The  computer  output  for  a sample  problem  using  TREAD/TLOOK  is  shown 
in  Appendix  C.  Both  of  the  previously  described  tables  are  input  but 
only  the  first  table  (table  reference  number  263)  is  used.  The  main 
program  used  to  implement  this  example  is  not  shown.  The  main  program 
only  calls  TREAD  to  input  the  data  and  TL00K  to  look-up  the  dependent 
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variable  values.  The  computer  output  resulting  when  TREAD  is  called 
from  the  main  program  is  shown  on  pages  C-2  through  C-5,  Appendix  C.  On 
pages  C-2  through  C-4  the  table  input  data  are  listed  as  processed  from 
the  input  cards.  On  page  C-5  a table  data  summary  is  shown.  The  purpose 
of  this  summary  is  to  show  the  user  how  many  tables  were  input,  what 
each  table  reference  number  is,  the  starting  array  location  of  each 
table,  and  table  storage  locations  used  and  locations  remaining.  Computer 
output  from  a main  program  using  repeated  TLOOK  calls  to  table  reference 
number  263  is  shown  on  pages  C-6  through  C-8.  These  data  are  plotted 
along  with  the  original  table  input  data  points  in  Figure  B-2.  As  shown 
in  Figure  B-2  the  interpolated  and  extrapolated  values  follow  the  curve 
very  satisfactorily. 
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APPENDIX  D 
TREAD/TLOOK 

DESCRIPTION  AND  FORTRAN  LISTING 


DESCRIPTION  OF  TREAD/TLOOK  CALCULATION  SEQUENCE 

The  input  section  of  TREAD/TLOOK  is  shown  on  cards  14  to  121, and  the 
interpolation  section  is  shown  on  cards  122  to  158,  Table  D-I.  In  the 
first  part  of  the  interpolation  section,  a search  is  performed  to  find 
the  correct  table  location  in  the  TDATA  array  to  begin  operating  on  the 
data.  On  cards  146  to  156  the  three  interpolation  or  extrapolations 
(one  for  each  independent  variable)  are  performed  using  the  SPLNQ1 
routine.  The  final  dependent  variable  value  returned  is  shown  on  card 
156. 
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DESCRIPTION  OF  FUNCTION  SPLNQ1  CALCULATION  SEQUENCE 

A brief  description  of  the  calculation  sequence  in  SPLNOl  is  as 
follows:  On  cards  8-20  of  Table  E-I,  the  storage  location  subscripts 

relative  to  the  beginning  location  of  the  input  data  in  the  X array  are 
computed.  On  cards  22  to  27,  the  given  independent  variable  value,  XIN, 
is  computed  with  the  low  and  high  value  of  the  input  independent  variable 
data  point  set.  If  XIN  is  greater  than  the  highest  data  point  or  lower 
than  the  lowest  data  point,  the  variable  NTRAP  is  set  to  0 or  1,  respectively. 
The  NTRAP  variable  is  used  as  a cue  for  extrapolation.  The  variable  L, 
on  card  30,  is  the  storage  location  of  the  upper  data  point  used  in  the 
previous  interpolation . L will  be  assigned  various  values  during  execution 
of  the  routine.  The  value  of  L is  initialized  zero  in  the  X array.  In 
the  situation  when  L is  equal  to  or  less  than  zero  a full  search  (cards 
34  to  39)  of  the  independent  variable  data  point  set  is  performed  in 
order  to  find  the  correct  interval  for  the  interpolation.  An  optimum 
search  technique  is  used  such  that  for  n data  points,  only  log2(n) 
comparisons  need  be  made  to  fir*''  he  interval  X(L  - 1)  XIN  <_  X(L) . 

In  the  case  of  large  arra  • . 50  points  for  example,  a maximum  of  6 
comparisons  need  be  made  to  find  the  interval  X(L  ~ 1)  f_  XIN  X(L) . 

In  the  second  situation  where  L is  greater  than  zero,  the  search  is 
begun  by  first  checking  the  data  point  interval  used  in  the  last  interpolation. 
After  the  correct  interval  is  found,  the  calculation  sequence  is  at 
cards  42  and  43.  The  variable  M which  satisfies  X(M  - 1)  XIN  X(M) 
is  stored  in  the  X array  for  the  next  interpolation  as  shown  on  card  43. 

On  card  46,  the  values  of  the  variables  L and  IQMODE  are  compared  with 
zero.  IOMODE  is  a user  set  variable  which  determines  whether  or  not  the 
second  derivatives  are  to  be  saved  and  used  in  the  next  interpolation 
sequence.  When  L = 0 these  data  points  are  being  interpolated  for  the 
first  time,  therefore,  the  second  derivates  must  be  calculated.  When 
IOMODE  = 0 these  data  points  are  being  interpolated  for  the  first  time, 
therefore,  the  second  derivates  must  be  calculated. 

The  condition  IOMODE  * 0 implies  that  in  each  usage  the  input  data 
points  to  SPLNQ1  are  changing  and  the  second  derivates  must  be  recalculated 
for  each  interpolation.  On  cards  47  to  86,  the  second  derivative  calculations 
are  performed.  The  second  derivative  calculations  are  optimized  to  the 
extent  that  if  IQMODE  = 0 only  the  second  derivatives  required  for  the 
current  interpolation  are  computed. 

The  calculations  on  cards  87  to  96,  concern  extrapolation.  A linear 
extrapolation  is  performed  using  the  slope  at  each  respective  end  point 
of  the  spline  curve.  The  linear  extrapolation  technique  is  preferred 
rather  than  a second  order  extrapolation  because  far  reaching  extrapolations 
are  more  controlled  and  predictable. 
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